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Abstract. We present a method of control which is able to create barriers to 
magnetic field line diffusion by a small modification of the magnetic perturbation. 
This method of control is based on a localized control of chaos in Hamiltonian systems. 
The aim is to modify the perturbation (of order e) locally by a small control term (of 
order e^) which creates invariant tori acting as barriers to diffusion for Hamiltonian 
systems with two degrees of freedom. The location of the invariant torus is enforced in 
the vicinity of the chosen target (at a distance of order e due to the angle dependence) . 
Given the importance of confinement in magnetic fusion devices, the method is applied 
to two examples with a loss of magnetic confinement. In the case of locked tearing 
modes, an invariant torus can be restored that aims at showing the current quench 
and therefore the generation of runaway electrons. In the second case, the method 
is applied to the control of stochastic boundaries allowing one to define a transport 
barrier within the stochastic boundary and therefore to monitor the volume of closed 
field lines. 
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1. Introduction 

Controlling chaotic transport is a key challenge in many branches of physics like for 
instance, in particle accelerators, free electron lasers or in magnetically confined fusion 
plasmas [llEllHlllliniEllZllHlinillDlinillSllH!. It covers a large variety of strategies. 
On the one hand, one can aim at recovering near integrable systems while on the other 
hand, one can request the existence of a boundary with strongly reduced transport 
that will act as some transport barrier. Another feature that will prove to be essential 
when implementing the control scheme is the relative cost of such a control device with 
respect to its merits. The cost that is considered here will be characterized by the 
magnitude and complexity of the Fourier spectrum of the control function with respect 
to that generating the spurious stochastic transport. In the case of systems such that 
it is essential to control a global transport property without significantly altering the 
original system under investigation nor its overall chaotic structure, one can restrict the 
control scheme to that of restoring a local transport barrier. In accordance with the 
latter idea, the control strategy that we develop here is based on building barriers by 
adding a small apt perturbation which is localized in phase space, hence confining all 
the trajectories. This local control strongly reduces the cost of the control scheme since 
only a small fraction of the phase space is modified. The counterpart is that the impact 
of such a localized reduction in transport must be efficient enough to be meaningful 
regarding the control strategy. 

In the case of particle accelerators, free electron lasers, and magnetic fusion, the 
magnetic field is designed to confine the charged particles. Hamiltonian dynamics allows 
one to describe the motion of the charged particles in such a magnetic field, and the 
confinement volumes readily take the form of invariant tori. The Hamiltonian formalism 
appears to be particularly well suited to investigate the control of such systems. In this 
article, we consider the class of Hamiltonian systems that can be written in the form 
H = Hq + eV i.e. an integrable Hamiltonian Hq (with action-angle variables) plus a 
perturbation eV. Provided the perturbation is not too small, the explicit dependence of 
eV on the angle variables breaks the invariance of the action variables and in a generic 
situation yields chaotic transport. The idea for controlling this transport is to modify 
the perturbation slightly and locally and create regular structures (like invariant tori). 

In the theory of magnetic confinement, the particle trajectories are characterized 
by three invariants, namely the magnetic moment, that is in fact an adiabatic invariant, 
the kinetic moment along the vertical axis which is the axis of symmetry of the 
axisymmetric equilibrium, and the particle energy (for a collisionless plasma). In steady 
state solutions, the various quantities of relevance for transport must solely depend on 
these invariants and not on the associated phases. When considering magnetic field lines, 
this six dimensional problem of the particle motion in phase space, can be reduced to 
a two dimensional one. The time variable then appears to be a toroidal angle, while 
the invariant is the toroidal fiux and the Hamiltonian, the poloidal fiux. The conjugate 
angle of the toroidal fiux is found to be a poloidal angle. In the case of a magnetic 
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equilibrium, the poloidal flux is solely a function of the action variable, the toroidal flux, 
both being labels of the magnetic surfaces, i.e. invariant tori. In this approach, and for 
an axisymmetric equilibrium, the energy, namely the poloidal flux, cannot depend on 
the angle since the latter depends linearly on the toroidal angle. It is therefore important 
that the unperturbed Hamiltonian for an axisymmetric equilibrium does not depend on 
the angle. 

In Refs. |141 IT^ ITHj . an explicit method of control was provided for generating / 
such that the controlled Hamiltonian Hq + eV + f is integrable. We point out that 
the control term we constructed is of order e^, i.e. much smaller than the perturbation 
which is of order e. This method of control has been applied experimentally to increase 
the kinetic coherence of an electron beam in a Traveling Wave Tube ^7j. One possible 
drawback of this approach is that the control term depends on all the variables and has 
to be applied on all the phase space. Here we provide a method to construct control 
terms / with a finite support in phase space, i.e. localized in phase space, such that the 
controlled Hamiltonian He = Ho + eV + e'^f has invariant tori whose explicit expressions 
are known. For Hamiltonian systems with two degrees of freedom, these invariant tori 
act as barriers to diffusion in phase space. 

In the present paper, this original method is applied to the control of chaotic 
transport in fusion plasmas. However, these examples are sufficiently general so that 
one can readily transpose the approach to other fields of physics. Let us consider two 
examples where the control of chaotic transport appears to be of particular interest. 

In a first example, one addresses the issue of large scale stochastic transport between 
two resonant surfaces. In tokamaks, this situation is met when two tearing modes 
develop and when their amplitude is large enough to satisfy the so-called overlap 
criterion proposed by Chirikov JH]- This leads to chaotic transport from one resonant 
surface to the other. The control scheme can target two results, first to maintain a 
closed magnetic surface (an invariant torus) so that confinement can be sustained until 
some other action is taken to suppress the tearing modes, second, should the tearing 
mode interaction have led to a pre- disruptive phase, to generate a transport barrier 
sufficiently robust to slow the current quench and provide conditions for a controlled 
disruption with no generation of runaway electrons [20] . 

In a second example, one considers the situation of an external modification of the 
magnetic equilibrium, for instance in the case of an ergodic divertor [21J where specific 
coils are implemented to generate chaotic transport in the outermost magnetic surfaces. 
A similar situation can be met in stellarators where the magnetic equilibrium can be 
such that a series of resonances in the plasma boundary creates a region of chaotic 
transport. In these cases, it can be of interest to control this chaotic boundary layer 
by restoring a closed magnetic surface at a given radius. In the latter case, one should 
consider a control scheme with similar properties to that generating the boundary layer 
with stochastic transport. This would ensure that a similar coil set to that generating 
the perturbation of the magnetic equilibrium can be used for the control purpose, hence 
ensuring some effectiveness in terms of cost and feasibility. 
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In Sec. 121 we provide an explicit construction of the localized control term. The 
equation of the created invariant torus is given explicitly. In Sec. El the formula of the 
control term for magnetic field lines is derived as well as the explicit formula for the 
magnetic surface which has been created by the localized control term. In Sec. |31 the 
two examples of magnetic perturbations are considered. It is shown on these examples 
that a small and apt modification of the poloidal flux is able to create a robust magnetic 
surface which prevents the diffusion of magnetic field lines. In the Appendix, we give a 
proof of existence of the control term. 

2. Control of chaos in Hamiltonian systems 

In this section we present a control method of Hamiltonian systems which is directly 
applicable to magnetic field lines. This version of localized control follows the one 
developed in Refs. [221I2S1- In order to understand the procedure to build the control 
scheme, let us first analyze the properties of a specific class of Hamiltonian systems that 
will provide the background for the control strategy. In the following A refers to the 
action variables and to the angle variables in a phase space of dimension 2L (hence 
both vectors are of dimension L). Let us consider Hamiltonian systems of the form 
H{A,6) = Ho{A) + eV{A,6). Using a suitable expansion, these Hamiltonian systems 
take the form 

H{A,e) = uj ■ A + ev{e) + w{A,e), (1) 

where {A,0) eW^ xT^,T being an angle space (torus) [— tt, 7r[ in a standard approach. 
The three contributions to H{A,6) are respectively Hq{A) = oj ■ A, the main term 
governing the integrable motion in the vicinity of A = 0, a perturbation ev{0) and a 
higher order term in A, w{A, 6) which can be written in the form 

w(A, 9) = ewi{e) ■ A + W2{A, 6), 

where W2 is quadratic in the actions, i.e. ^2(0, ^) = and 9aU'2(0, 0) = 0. 

The vector u: contains the frequencies of the quasi-periodic motion defined by A = 
and 6 = uj. Furthermore, a; is a non-resonant vector of M"^, i.e. there is no non-zero 
k G Z"^ such that ■ k = 0. Without restricting the computation of the control term, 
one can assume that w{0, 0) = for all eT^ and that /^^ v{6)d^6 = (the average 
of V is set to zero). We consider a region near A = 0. We notice that for e = and for 
any w not necessarily small, the Hamiltonian H has an invariant torus at A = 0. The 
controlled Hamiltonian we construct is given by 

ife(A, 6)=u:-A + ev{0) + w{A, 6) + E^f{e)VL{A, 0), (2) 

where is a smooth characteristic function of a region around a targeted invariant torus. 
We wish to stress that the actual control term e'^f that we construct only depends on 
the angle variables 6 and is of order e"^ (see Appendix) . We prove that the control term 
is : 

f{e) = -e-M-eTdev,e), (3) 



Control of stochasticity in magnetic field lines 



5 



where dev denotes the first derivatives of v with respect to 6 : dev = X^kez^ i\ivi^e^^'^ 
and where the linear operator F is the inverse operator of the convective shift 6 ■ de 
relative to the unperturbed dynamics governed by Hq{A). Its explicit form acting on a 
function v{0) = Xlkez^ Vy^e^^^'^ is : 

Yv{0)= y -^e'^', 

keZL\{0} 

The operators T and de commute so that d0{Tv) = T^dev). 

There is a significantly large freedom in choosing the function Q. It is sufficient 
to have Q{A,0) = 1 for ||A|| < e. For instance, Q{A,0) = 1 would be a possible and 
simpler choice, however representing a long-range control since the control term f{6) 
would be applied on all phase space. On the opposite way, we can design a function 
Q such that the control is localized around the created invariant torus : We denote 7^ 
and Tp two neighborhoods (in phase space) of the targeted invariant torus such that 
Ta d Tj3. The characteristic function is chosen such that = 1 if x G 7^ and 

Vt{x) = if X ^ 7g, and VL is smooth on all phase space. We choose the characteristic 
function Vt to be 

n{A,e) = VLU\\A + eVdev\\), 

where VLioc{x) = 1 if x < a and f2/oc(x) = if x > /9 and, e.g., a polynomial or another 
function for x e]a, /3[ such that VLioc is smooth on M+. The main advantage of this step 
function Vt is that the control needs less energy (only in the part of phase space where 
the control is localized) and also it does not change the other part of phase space. 

In these cases, we prove that He given by Eq. Q has an invariant torus located at 
A = —eVdev. For Hamiltonian systems with two degrees of freedom, such an invariant 
torus acts as a barrier to diffusion. For the construction of the control term, we notice 
that we do not require that the quadratic part of w is small in order to have a control 
term of order e^. Moreover, if the lowest order in powers of A of w is n > 2 then the 
control term is of order e". 

In order to derive the expression of the control term, we consider the canonical 
transformation generated by 

F(A',6>) = e-A'-eTv{0), 

which is a translation in action by eVdev, i.e. {A',0') = (A + eVdev^O). The 
Hamiltonian is mapped onto 

i^,(A', e') = i^-A' + w{A! - eVdev, 0') + ^^(^O^dl A'||). 

The translation in action is such that the contribution ev is canceled out since u) ■ Tdev = 
V. The action A' = is found to be conserved by the flow of He since 

dA' 

= -dew{A' - eVdev, 6') + e{Vdlv)d^w{A' - eVdev, 6') - e\dem{\W\). 
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and since the control term f{0) defined in Eq. Q, is such that : 

def = -e-^dewi-eVdgv, 6) + e-\Tdlv)dAw{-eTd0V, 6). 

As a consequence, in the domain such that f2(||A'||) = 1, we find A' = imphes 
dA'/dt = 0, and consequently A' = is an invariant surface for He- 

3. Magnetic field lines 

The magnetic field line dynamics in a toroidal geometry can be written in a Hamiltonian 
form [211 EHl El 

de _ dH 

dip dip ' 
dijj _ _dH 

~d^ ~ ~~de' 

where ip which plays the role of effective time is the toroidal angle, if) is the normalized 
toroidal fiux and H is the poloidal fiux. The poloidal angle 9 is the conjugate variable 
to the action ip. We consider the following class of Hamiltonian systems : 

6, if) = HoW + sH.ii:, 9, if). (4) 

We denote Q{ip) = Hq{iP). The quantity q = 1/Q is the safety factor. For e = 0, 
we recover the unperturbed magnetic equilibrium such that dip/dip = 0, ip = ipo are 
invariant tori also characterized by the rotational transform q{ipo) = dip/dO. We select 
a magnetic surface by its unperturbed action ip = ipQ where one wants to build a barrier 
to diffusion. We expand Hq into : 

CO _ 

We denote u = Q{ipo) = l/'?('^o) the winding ratio of the selected magnetic surface, and 
Q*^'^ denotes the l-th. derivative of Q. Following the notation of Sec. |2l we have 

v{e,(p) = Hi{ipo,e,ip), 

and 

wiip, e, if) = HoiiP) - Qm{ip - ^o) + e\E^{iP, ip) - H.iiPo, 9, y.)]. 
We expand w : 

oo ^ 

w{ip,e,v) = ed^H,{iPo,e,^){iP-iPo)+J2-^^—^_[Q^'\^^^ {iP-iPoY 

We notice that w{ipo,0,ip) = and that d^w{ipo,9,ip) = ed^Hi{ipo,9,ip) is of order e. 
Thus this Hamiltonian satisfies the requirements to construct a localized control term 
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Following Refs. j2Zl 123 123 , the Fourier expansion of Hi is given by 



Hi{ipo, 6*, V?) = ^ Hmnii'o) cos{m9 - rup + 



for some constant phases Xmn- In order to apply the control procedure described in the 
previous section, we consider that ip is an angle variable and E its conjugate action. 
In this way, we map the "non-autonomous" Hamiltonian system with 1.5 degrees of 
freedom into an "autonomous" Hamiltonian of the form with two degrees of freedom 
where the actions are A = {ip — ipQ, E) and the angles are 6 = [6, if). Hamiltonian P|) 
has the form (0) with u: = (tu, 1). 
The control term is given by 

f{e, = [d^HiitPo, 0, V)] TdeHiiiJo, 0, ^) 

- E TTXTm (^^'H^o) + ed'^+'Hii^o, 0, ^)) {TdeHii^o, 0, y.))'^^ , (5) 
1=1 ^ + ^■ 

where 

TdeHi{ipo, e,ip) = y^ mHmn{i^o) ^^^^^q - mp + Xmn)- (6) 

•i— ^ mi,] — r) 



m,n 



Using this control term, the controlled Hamiltonian has the invariant torus whose 
equation is 

^ = ^o-eTdeHi{iJo,e,ip), (7) 

The difference between and ipo, of order e, is a function of ipo both via the dependence 
of Hi on ipo and that of the operator F with respect to the frequency u = Q^ipo). The 
dominant term of the control is given by considering only / = 1 in Eq. (jSj) , which leads 
to 

f2{e,v) = [d^Hi{iJo,9,^)]TdeHi{^o,0,^)-^^^{TdeHi{iJo,0,v^^^ 
where TdeHi is given by Eq. (jHl). The full control term is given as a series in e : 

oo 
s=2 

where fs is given by 

fs{e,^) = p%id;-'Hi)iTdeHir-' + ^zll!lig(-l)(^,)(^9,/70^ (9) 
[s — 1)1 ^ si 

for s > 2. 

In summary, the controlled Hamiltonian we consider is 
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where / is given by Eq. for the complete control term, or by Eq. (jSI) if one wants to 
use only the dominant term, and TdgHi is given by Eq. ®. The result in Sec. El ensures 
that the controlled Hamiltonian He with the control term / given by Eq. ^ has an 
invariant surface whose equation is given by Eq. ((7j). These expressions provide the basis 
to investigate the control of the stochastic transport governed by magnetic perturbations 
in fusion devices in the two cases, core tearing modes and boundary stochastic layers as 
done in Section 4. 



4. Control of stochastic transport in fusion devices 

4-1. Control of the loss of confinement governed by coupled tearing modes 

We consider a Hamiltonian system where ip is the normalized toroidal flux, 1^ = 1 
at the plasma boundary, with a g-profile as chosen in Refs. [SUl EH 121] 

In this expression, q{ijj) is a monotonic function of ip that varies between 1 on axis, for 
ip = 0, and 4 at the plasma edge, -0 = 1. With such an expression, the slope d\og{q) / dip 
is approximately constant and equal to log(g('?/' = 1)). The results do not depend on 
the precise form of the safety factor profile and can be readily extended to any profile. 
We are interested here in the case of two tearing modes with low rational numbers, 
in practice (m, n) = (3,2) and {m,n) = (2,1) where m is the poloidal mode number 
and n the toroidal mode number, locked to one another. In the case analyzed here 
both perturbations have equal magnitude characterized by e. The latter parameter is 
assumed to be small, but such that the onset to stochastic transport is reached with a 
magnetic perturbation of the form 

Hi{ip, 9, ip) = cos(2e + cos(30 - 2ip), (12) 

so that the Hamiltonian of the system is 

^ = / -Tk + eH,iip,e,^). (13) 

The resonant surfaces are found to be located on iIj^ 2 ~ 0.266 {q{jpz,2) = 3/2) and 
V'2,1 ~ 0.456 (g(V'2,i) = 2). Expanding to second order the Hamiltonian H in the vicinity 
of the two resonant surfaces and retaining the resonant term of ifi, one recovers the 
characteristic Hamiltonian of the pendulum [S2] that allows one to define the so-called 
unperturbed island width 5 in units of ip : 



5 = 2 



1/2 

e 



-d{l/qm 



A Poincare section of the dynamics given by Eqs. (jl3j) - (jl2j) is represented on Fig. Qfor 
e = 0.004. For this value of e, the island widths are 53^2 ~ 0.125 and ^2,1 ~ 0.147 so 
that a Chirikov parameter a = {62^1 + ^3,2)/(V^2,i — ^^3,2) can be computed and is about 
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>Q.35 




Figure 1. Poincare sections of Hamiltonian H13(l witli Hi given by Eq. H12|l with 
e = 0.004 for V' e [0,0.7] (left panel) and for V e [0.2,0.5] (right panel). 



1.4, hence larger than the reference value for the onset of large scale transport between 
the resonant surfaces '?/'3^2 and iIj2,i as readily observed on Fig. ^ However, e is small 
enough that significant remnant islands are still present. 
The control term we apply is given by Eq. (jHJ where 

d^Hi{ilJo,0,ip) = 0, 

2 . 3 



TdgH,{^o,e,ip) 



cos(26' -(f) + 



cos(30 - 2ip). 



2uj-l ' " 3cu - 2 
For t/'O) we choose ^/'o = 0.35 hence between the two resonant surfaces ^/'3_2 ~ 0.266 and 
^2,1 ~ 0.456. Choosing other values of is equivalent to moving the barrier one wants 
to create. The expression of the partial control term /2 is given by Eq. (jH)) : 



2 cos(2^ -v) ^3, cos(3^ - 2ip) 



(14) 



■tp=ipo V 2uj — 1 30; — 2 

where uj = 1/ q{ipo) = (2 — ^/'o) (2 — ^ipo + t/^q) /4 and where we used the fact that Hi does 
not depend on ip. For the present value of ipo one finds g(V'o) ~ 1-7 and thus u ~ 0.587. 
The full control term / creates an invariant torus whose location is given by 

2 , 3 



ip = ipo 



cos(20 — ip) 



cos(30 - 2ip) 



(15) 



^2u;-l ' 3uj-2 
The magnitude of the angle modulation of the invariant torus labeled by ijj is relatively 
large (larger than 40e). 

For the localization function Q, we use two choices : The first one is f^(|x|) = 1 
for all |a;| G M^, i.e. without localization. Such a control procedure appears to be more 
readily applicable in fusion plasmas where it is difficult to device an electromagnetic 
perturbation localized in a narrow region of the core plasma. The Poincare sections of 
trajectories of the controlled Hamiltonian Hq + eHi + are represented in Fig. El for 
trajectories started from below or from above the invariant torus given by Eq. ()15p and 
represented by the bold curve. We notice in particular that the trajectories started from 



Control of stochasticity in magnetic field lines 



10 




Figure 2. Poincare sections of Hamiltonian 1)13(1 with the control term /2 given by 
Eq. p4(l with e — 0.004 and i/jq = 0.35 using il = 1, tlius witli no locaHzation of 
the control term : with initial conditions below (left panel) or above (right panel) the 
surface given by Eq. 1(15(1 plotted in bold. 



below (resp. above) the invariant torus remain below (resp. above) it. The impact of 
the control term is noticeable away from the target torus ipo since one can observe much 
larger remnant islands in the vicinity of the two resonant surfaces compared to the case 
without control. 

The second choice of f2 is a localization function. Although such a localized control 
term has yet no clear application to fusion plasmas, it demonstrates that the efficiency of 
the control procedure is not tied to the perturbation of the whole phase portrait but only 
to a specific region of the phase space defined by Eq. (fT3|) . It is therefore illustrative to 
consider such a control procedure in order to check in particular that the leading effect 
of the control scheme is restoring a selected magnetic surface and not a wide range 
cancellation of the effects of the perturbation. We choose flioc{\Sip\) = 1 for l^i/^l < 6ipa, 
^Zoc(I^V'l) = for \Sip\ > Sipp and a third order polynomial for \6ip\ ElSipa, Sil'/si for which 
Qioc isaCi-function,i.e. niUM) = l-im-S^|Jafi^S^Pfs-6^^-2\6^P\)/i6^|Jf,-6^|J^)^ 
In principle, one can choose arbitrarily small values for Sipa and 6ipi3 if one uses the full 
control term / given by Eq. (jSJ. However, since / is given by a series, it is more 
practical to consider the truncated control term /2 (or a truncation of the series which 
gives the control term /). Then the value of 6tpa has to be not too small such that the 
set {{A,0) s.t.f2ioc(A, 0) = 1} contains the invariant torus which, for /2, is e^-close to 
the one obtained using the complete control term. This leads to a restriction for Sijja, 
5ipa ^ 10~^ to embed the lowest order in the control term expansion. For the numerics, 
we choose = 0.01 and 5il)p = 0.02. A Poincare section of the dynamics of He given 
by Eq. (jlOj) with the control term /2 is represented on Fig. El for e = 0.004 and for 
Q = Qioc- The bold curve corresponds to the invariant torus given by Eq. (fTK|l that has 
been created by the addition of the control term which is localized around this surface. 
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From these figures, we clearly see that the upper and lower parts of phase space are 
very similar to the ones of Fig. ^ (without control). More precisely, we notice that the 
structure of the resonant islands is not modified, even the neighboring ones. What has 
changed is the dynamics in the neighborhood of the bold curve due to the action of the 
localized control. There is now an invariant torus which prevents the motion to diffuse 
from the lower part to the upper part of phase space. These two parts are invariant by 
the dynamics of the controlled Hamiltonian. 

The invariant torus created by the localized control persists to higher values of e. For 
e > 0.1, the trajectories start to diffuse through the broken invariant surface. The 
diffusion can be reduced or even suppressed by taking into account higher order terms 
fs for s > 3 in the control term series. We point out that the value of e for which the 
partial control term /2 is efficient depends on the choice of ipQ. There is freedom to 
choose the initial surface ipo provided that q{ipo) is sufficiently irrational. 

Let us define the norm of a function f{9,ip) = Ylimn fmnCos{mO — mp) as ||/|| = 
maXm,ra \fmn\- If wc compare the relative size of the control we obtain : 

typically 14% of eHi for e = 0.004. The magnitude of the control term computed in 
such a way appears to be a small fraction of the magnitude of the initial perturbation 
that led to the stochastic transport. Moreover, when considering the localized control 
term, the control only acts on a finite and small portion of the phase space [0, 1] x 
around the invariant surface (the size of the support of flioc is 4%). 

Among the aims of a control scheme, that of restoring a region with closed magnetic 
surfaces in a disruptive phase can prove to be very valuable even if all nested magnetic 
surfaces are not recovered. Indeed after the thermal quench, when most of the plasma 



Control of stochasticity in magnetic field lines 



12 



kinetic energy has been lost, the standard disruption scenario enters the so-called current 
quench. In this phase the plasma current strongly decreases which leads to a loop- voltage 
spike in order to sustain the plasma current. Such a voltage spike creates conditions 
for the transfer of the energy from the poloidal system to the plasma and specifically 
for the generation of runaway electrons. These are considered as a major problem in a 
device like ITER since these runaway electrons will lead to very strong energy deposition 
on reduced areas and deep into the material given the slowing down distance of these 
energetic particles. As a consequence, this phase of the disruption can lead to serious 
damage to wall components. Defining appropriate control tools to reduce the energy 
transfer to a runaway population or to lower the runaway population is thus a key 
task for ITER. One can readily assume that already in the thermal quench phase, a 
broad spectrum of magnetic perturbations is responsible for the confinement loss due 
to the onset of large scale transport governed by the overlap of the various modes of 
the spectrum. Among these modes, one can expect that the low (m, n) modes, such as 
the (3,2) and (2,1), to overlap and govern a strongly enhanced transport as required 
for the thermal quench. When applying the control magnetic perturbation towards 
before the onset of the current quench, very likely when the stochastic core region 
comes into contact with the walls, one aims at restoring a magnetic surface between 
the two resonant surfaces. In practice, this yields a transport barrier (with respect 
to the very degraded confinement during the thermal quench) that will separate the 
core region, where the total current will be maintained, from the edge region where 
all the current will be lost due to the connection to the wall components. Such an 
insulation ought to prevent the current quench and the generation of runaway electrons. 
Should the loop-voltage spike be only partially suppressed, and a fraction of runaway 
electrons still generated, these should remain confined within the transport barrier. The 
control scheme applied to a tokamak disruptive phase thus appear to provide conditions 
for a soft landing with no current quench. Applicability of such a scheme depends of 
course on the possibility to generate the spectrum of the control term. Let us consider 
the spectrum of /2 as given by Eq. (fT^ . It consists of four modes, the two modes 
corresponding to the tearing modes that define the region where a magnetic surface is 
to be restored, hence the (3, 2) and (2, 1) modes and their nonlinear combination, namely 
the (1, 1) and (5, 3) modes. The relative magnitude of the modes determines the value 
of uj and therefore the location of the barrier in terms of the safety factor q. For a given 
discharge, this parameter can be preset so that the activation of the control scheme will 
not require a real time computation of the required spectrum. The means to generate 
such a magnetic perturbation is to implement dedicated coils at the plasma boundary. 
Of course, this means that one must contemplate the non-localized control. Although 
it is not the purpose of the present paper to design the appropriate coils to generate 
these modes, one can underline the main features of the coil system. First, the required 
coils are significantly simpler than the coils required for boundary plasma control either 
of the ergodic divertor kind [SHI Olj or to control ELMs ^I] . This is due to the low 
poloidal mode number of the modes that will therefore exhibit a moderate radial decay. 
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Furthermore, the magnitude of the control term is typically 14% lower than required 
to generate stochasticity as required for boundary plasma control. Of course a proper 
project would have to address many issues left open here, like the effect of a mismatch 
in magnitude of the control term, or the effect of spurious modes on the overall control 
efficiency. However, the tests that have been performed to examine the robustness of 
the control procedure (see Refs. [T3J [T3]) give us confidence that these effects should 
only lead to moderate reduction of the control efficiency. These favourable features 
show that the application of the control scheme to prevent the current quench during a 
disruption should not face outstanding design difficulties especially when compared to 
the consequences of uncontrolled runaway generation. 

4-2. Control of stochastic boundary plasmas 

For the second example, we consider a magnetic perturbation which models the magnetic 
field lines in an ergodic divertor [HllEIlEHllSSlEn]- We use the g-profile given by Eq. (|TT|l 
and the following magnetic perturbation 



where n = 2, 6^ = vr/3, mo = 6 and the sum ranges from mo — 4 to mo + 4. We 
have chosen m = m. In experiments, the strong modulation of the toroidal magnetic 
field between the inner region and the outer region of the torus can lead to a significant 
departure from this approximation. This is the case for the Tore Supra experiment with 
m ~ m/2 and the Textor DED experiment with m ~ 2m However, such aspects 
of the perturbation, while important for the engineering constraints, do not modify 
significantly the computation of the control term. Also, the values chosen for the mode 
numbers are relevant for the Textor DED, a rather low m and n configuration that is 
required to provide a reasonable stochastic boundary in spite of the DED coil location on 
the high field side jHSj- A possible operating regime of the DED is to induce a rotation of 
the magnetic perturbation at constant frequency for all the modes of the spectrum by an 
appropriate phasing of the coils. The control scheme derived here does not depend on the 
choice of the rotation frequency provided the magnetic perturbation of the control term 
is designed to achieve the same rotation frequency as the main magnetic perturbation. 
For the purpose of the present analysis, the rotation frequency is not taken into account, 
and we shall concentrate on the control of a stochastic boundary induced by a steady 
state perturbation. 

Few modes are then effective in the spectrum given by Eq. ()lfi|l and the largest 
extent of the perturbed edge region is achieved when the main mode of the spectrum, 
mo = 6 and n = 2 (hence q{ipG,2) = 3) is located within the edge region, at ?/'6,2 ~ 0.75 
for qedge = 4 for the chosen profile of the safety factor. The region chosen for the control, 
V'o ~ 0.9, q{ipo) ~ 3.6 thus appears to be located in the vicinity of the mode (7, 2) of the 
spectrum characterized by a large amplitude (~ 0.8 to be compared to 1 for the main 
component (6, 2)). 




m 



(16) 
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Figure 4. Poincare sections of H given by Eq. witli given by Eq. 1)16(1 
e = 0.003 for V' e [0.5, 1.1] (left panel) and for V e [0.71.1] (right panel). 



with 



One well documented operational limit of stochastic boundary plasmas is achieved 
when the stochastic region extends over a too large radial extent jSHlEZllSHI and when the 
stochastic boundary actually reaches the wall, namely when all invariant tori between 
the stochastic region and the wall are destroyed jHEl EH EH] • A Poincare section of the 
dynamics is represented on Fig.|3]for e = 0.003 where we notice that magnetic field lines 
diffuse up to outer-edge of the plasma {ip = 1) and to the wall. 

The control term we apply is given by Eq. (jSJ where 

A tr M a \ sm[{m-mo)ed]m ^/2-i , . 

u f I a \ .^„^sin[(m-mo)^d] m/2 m 
TdeHi{%l>Q, e,(p) = } (-1) — ^pQ cos(m6' - nip). 

m ^ ' 

The control term thus appears as the product of two terms (given above) each having 
a similar spectrum as the initial perturbation given in Eq. As a consequence of 

this nonlinearity all coupled modes m\ + m2 and n\ + n2 must be generated leading 
in practice to a very large spectrum for the control perturbation. Since the original 
spectrum is in fact induced by a modulated magnetic perturbation with mode number 
mo, no localized in a finite poloidal arc (that govern the width of the spectrum), the 
control term will, at first order, have to be produced within the same poloidal arc but 
with a modulation 2mo, In^- The required coil is therefore very similar but with larger 
main mode numbers. The design of such a coil in the case of the DED would not 
be possible due to specific constraint of this design while it would have been possible, 
although difficult, to implement the required perturbation with the Tore Supra ergodic 
divertor [HHl El]- In the case of the DED, it is thus of interest to consider the situation 
of a degraded control term such that only a subset of the optimum control term is 
implemented. Before analyzing this aspect of the control procedure and its impact on 
the coil design, let us consider the effect of the full control term and its impact on the 
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Figure 5. Poincare sections witli Hi given by Eq. Hlfi|l using — ilioc, with the control 
term 2^/2 given by Eq. JHl with e = 0.003 and ipo ^ 0.92 : with initial conditions below 
(left panel) or above (right panel) the surface given by Eq. ijjj). 



transport within a stochastic boundary. 

We use the same characteristic function Qioc as in the first example with Sipa = 0.02 
and Sipp = 0.03. We choose ^po ~ 0.92. A Poincare section of Hamiltonian (fT^ with Hi 
given by Eq. (fTT)|l and with the control term /2 given by Eq. (jH)) for this example with 
e = 0.003 is represented on Fig. El 

As for the first example, we clearly see that the upper and lower parts of phase 
space are very similar to the ones of Fig. |3] In other words, this localized control does 
not affect the diffusivity of magnetic lines above the magnetic surface around ip = ipQ. 
However the two parts are disconnected by the dynamics since they are invariant by the 
dynamics of the controlled Hamiltonian. A magnetic surface whose equation is given 
by Eq. ((Zj) has been created and acts as a barrier to the diffusion toward the border 
of the plasma. As stated above, there is experimental evidence that indicates that 
such a transport barrier is sufficient to decouple the core plasma from the edge plasma 

[SSI EH EH]. 

The norm (as defined in the previous section) of the control /2 is about 15% of eHi 
for e = 0.003. Moreover, the control only acts on a finite and small portion (around 
7%) of the phase space [0, 1] x around the invariant surface. 

For 6 > 0.003, the truncated control term /2 is not sufficient to create the barrier. 
Therefore, one has to take into account more terms in the series expansion of the control 
term. A Poincare section of Hamiltonian pOj) with Hi given by Eq. p6p. with the 
control term and with the control term e'^f2 + e^fs given by Eq. for e = 0.004 
is represented on Fig. IHl 

We clearly see that /2 is not sufficient to create an absolute barrier to diffusion in 
contrast with the control term obtained with the addition of /s. For higher values of 
e, more terms of the series can be included if necessary. However we notice that there 
is still an effective barrier to diffusion which prevents most of the trajectories to diffuse 
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Figure 6. Poincare sections with Hi given by Eq. (|16|l with the control term 2^/2 (left 
panel) and with + e^fs (right panel) given by Eq. with e = 0.004, V'o ~ 0.92 
and VL = \, and with initial conditions below the surface given by Eq. l(7|l. 



0.015 




Figure 7. Probability distribution functions of trajectories in the presence of the 
partial control /2 (dotted line) and the partial control /2 + /s (bold line). 



toward the edge. In order to measure the efficiency of the control in the case where the 
partial control term /2 is not sufficient to create a barrier, we have plotted in Fig. [3 the 
probability distribution function of trajectories launched from below the barrier (|7j) in 
the chaotic sea for a fixed elapse of time. This function is averaged over the angles d 
so that the barrier is not strictly localized on ifjQ = 0.92 due to the angle dependence 
(set by Eq. ((Zj)). We notice that /2 is still efficient to reduce the diffusion of trajectories 
since most of the trajectories are such that ip < 1. 

Moreover, there is the freedom in simplifying the control term by removing some 
Fourier components. This can be analyzed when considering the control term computed 
for the interaction of two tearing modes (see Eq. ()14|) in Sec. 14.11) . and generalized to 
any set of modes, {mi,ni) and (m2,n2) with amplitude Ai and A2 characterized by a 
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weak dependence. In this case, the control term is roughly : 



X 



tp=t(!0 

/ rriiAi cos{mi9 — riiip) ^ 7712^2 cos(m26' — n2f) 
\ miuj — Til m2L0 — n2 

Such a control term is characterized by four Fourier components, two are corrections 
to the amplitude of the original components and two are coupling terms with angle 
dependence (mi + 7712)6* — (ni + n2)v? and {mi — 1712)0 — {rii — n2)(p. Only the first mode 
is resonant between the resonant surfaces of main components and has therefore the 
largest weight in the control procedure restoring an invariant torus in its vicinity. This 
analysis can be transposed to the control of the stochastic boundary addressed here, 
hence such that ni = n2 = n. For instance if one wants to create a magnetic surface 
between the island with period 7 and the one with period 8, the Fourier mode with wave 
vector (15,4) is dominant. In the case of a localization between resonant islands with 
frequency vector {mi,n) and {m2,n), the approximate control term then reduces to 

f2,liO, if) = frmm2 COs[(?7li + 77^2)6' - 2nLp], (17) 



with 



f _ / 1 ^m^+m2 sm[(mi - mo)grf] sin[(m2 - mo)grf] ^^(^,+^,)/2 

Jmim2 \ ^) I \ I \ ro 

7r(?77i — mo) vr(m2 — mg) 
mim2 / (mi + m2)ijJ — 2n <5'(V'o 



X 



{rriiuj — n){m2UJ — n) \ Aipo 2 

In this expression, the ip dependence of the perturbation Hi is taken into account leading 
to a correction in the magnitude of the control term. Of course, the theorem does no 
longer ensure the existence of the invariant torus since this control term is approximate. 
However this simplified control term requires less energy (the ratio between the energy 
necessary for the control and the one of the magnetic perturbation is 4% for e = 0.003). 
Furthermore, the strong simplification of the spectrum of the control term should 
translate into the design of the dedicated control coil. The effect of the simplified 
control term of Hamiltonian (jl6p given by Eq. (jl7|) can be seen on the Poincare section 
of Fig. ISl It clearly shows that an invariant torus bounding the motion of magnetic 
field lines has been created by this simple control term for e = 0.002. However, when 
e is increased to e = 0.003 there is no longer an invariant torus and field lines leak 
out toward the ^/^ ~ 1 values. The density of points in the Poincare section of Fig. |H1 
is indicative of the existence of a transport barrier that inhibits the transport at the 
location of the invariant torus observed when the full control term is applied. 

In a similar way as in Fig. IHl we have measured the efficiency of the control in the 
case of Fig. IHl when e = 0.003 by plotting in Fig. El the probability distribution function 
(PDF) of trajectories launched from below the barrier (|7j) for a fixed time. The bold 
line represents the PDF with the control term /2 which is sufficient to create a barrier 
to diffusion when e = 0.003. The thin line is the PDF without control and the dotted 
line represents the PDF with the simple control given by Eq. (fTTjl . We notice that the 
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Figure 9. Probability distribution functions of trajectories without control term, 
with the simplified control term (|17|) (dotted line) and with the control term /2 (bold 
line). The dashed lines indicate the V-'-extcnsion of the created invariant torus. Here 
e = 0.003. 



simplified control term is still efficient to reduce the diffusion of trajectories (by a factor 
2 for the value of the PDF). In this strongly simplified control term is found to 

provide the required reduction of transport on the prescribed surface. 

Regarding transport issues in a fusion device this remnant loss of magnetic 
confinement can be comparable to the existence of an invariant torus provided the 
transport across the barrier along the field lines is smaller than the turbulent plasma 
transport jS^ . 

In this Section, we have shown that a transport barrier consisting of an invariant 
torus or a region of reduced field line diffusion can be generated within a stochastic 
boundary. This possibility can prove to be important not only for ergodic divertor 
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experiments as addressed here but also for other fusion device such as stellarators, 
where such stochastic boundaries are intrinsic to the magnetic equihbrium for a set 
of configurations jlOl 1^ as well as Reverse Field Pinches where restoring the 
magnetic surfaces in the outermost region is a key to enhanced performance jl^ . Finally, 
the present results can also be used to analyze the ELM control scheme addressed on 
DIII-D jUj. In the latter experiment a very weak magnetic perturbation is shown to 
control the large ELMs associated with the H-mode operation. Analysis of the magnetic 
structure have indicated that the magnitude of the perturbation is comparable to the 
error field due to toroidal coil misalignment. There is therefore a possibility that a 
control of stochastic transport of the type analyzed here is at work in these experiments. 
To investigate such a process one would need to reach a precise modeling of the ELMs, 
as proposed in Ref. to investigate the impact of such a control on the ELMs and 
H-mode pedestal physics. This is clearly beyond the scope of this paper. 

5. Conclusion and discussion 

We have developed a strategy of localized control of Hamiltonian systems. This control 
scheme is original insofar that rather than compensating for the perturbation driving 
the stochastic transport it is based on adding an extra perturbation of lower magnitude 
with respect to the initial perturbation in order to restore an invariant torus. A simple 
demonstration of this control scheme is proposed in Sec. El The application to magnetic 
field lines, in Sec. El allows one to show that it is possible to create isolated magnetic 
surfaces within a region of stochastic field lines. This prevents magnetic field lines 
from diffusing throughout the stochastic region and thus generates two independent 
regions separated by a transport barrier. It is important to stress that in such cases 
the perturbation leading to stochastic transport is of order e while the control term 
is of order e^. We note that the control we propose can be applied in all generality 
since the construction is independent of the specific form of the Hamiltonian. Moreover, 
this method is also applicable if the spectrum of the magnetic perturbation is given 
numerically. Control for two cases relevant to plasma confinement in magnetic fusion 
devices are addressed in Sec. |31 A rather generic case of the control of stochastic 
transport generated by two neighboring magnetic islands is provided by the overlap of 
low mode number tearing modes. In this case, the control term is readily computed and 
provides a control scheme both in the case where the control term is localized on the 
area of the restored magnetic surface and in the more relevant case of a control term 
acting on the whole plasma. The latter control scheme is more in line with a control 
perturbation generated by external coils. Such a control can prove to be important to 
confine the plasma current during the current quench phase that is readily associated 
with stochastic transport governed by coupled tearing modes. The second example 
is the control of stochastic boundary layers. In such a situation, it can prove to be 
important to control the region in contact with the wall via the stochastic boundary. 
The calculation of such a control scheme is performed in the ergodic divertor framework 
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and can readily be extended to other configurations such as stellarators. In this case, 
it is shown that a Fourier truncation of the control term still provides a control of 
transport at the chosen location in the boundary layer. This underlines the robustness 
of such a control scheme. It is interesting to note that the control term corresponds 
to the second harmonics of the spectrum used to generate magnetic perturbations. In 
a standard analysis, these harmonics are not taken into account. There is therefore a 
possibility that those parts of the spectrum that are usually neglected, because of their 
low amplitude and higher harmonics play a role in reducing locally the transport across 
the stochastic boundary. Going one step further, it is possible, using the ideas presented 
in this paper, that a weak magnetic perturbation generated by external coils acts as a 
control term of a boundary magnetic perturbation due to coil misalignment. This issue 
was raised after the first tests of ELM control on DIII-D with the so-called I coils |4nj . 
The present results indicate that such a counterintuitive property cannot be ruled out 
without considering the possibility of checking that the controlled perturbation and 
intrinsic perturbation are not related according to the relationship presented in Eq. (jHl). 
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Appendix A. Existence and regularity of the control term 

In this Appendix, we prove that the control term and the canonical transformation are 
bounded. It reduces to prove that Tdgv is bounded. This follows from usual KAM 
proofs (see for instance Ref. [IH]) which also shows the existence of invariant tori. The 
Fourier expansion of Tdev writes : 

^-^ ijj ■ k 

kezi\{o} 

We assume that a; satisfies a Diophantine condition, i.e. there exist a > and r > L — 1 
such that 

leu ■ kp^ < (T||kf , VkGZ-^\{0}, 

where ||k|| = \k^. Moreover, we assume that v is of class (7''+'^+^ where r > 1, i.e. 
bounded for the norm of scalar functions v{&) = Xlkez^ fke*'''^ : 

= |t;o| + ^ |t;k|||k||''+^+\ 
kez^-Xfo} 

In the same way we define the norm of vectorial functions of class as ||g||r = 
Z^ili \\9i\\r- Thus we have 

\\rdev\\r= 

kezi\{o} 



llkll 


l^kl 




•k| 



Control of stochasticity in magnetic field lines 



21 



The Diophantine condition of lo gives 

\\Tdev\\r Yl l|kr+'+1'^^k| < (T\\v\\r+r+l. 

This norm is bounded and hence Tdev is of class . Hence the loss of regularity 
between v and Tdev is the constant r + 1. We notice that we could have also weakened 
the hypothesis on {v,u) into Hr^efUr < oo. 

In order to have an estimate on the control term, we denote 

r]i = sup WTdevW < a\\v\\r+r+i, 

r]2 = - sup sup max \dAiW{A,6)\. 

If w is given by w(A, 0) = ewi{6) ■ A + 1^2 (A, 6), we see that 

^72 < ||wi||r + r/i sup sup ma.x\dA,dAj'W2\. 
\\A\\<sm eeJL i,j 

The control term given by Eq. © can be rewritten as 

= / dA.w{—seTd0V,6) ■ eVdevds. 
Jo 

Thus we have 

sup 1/(0)1 < r]iV2, 

where rji and 772 are of order 1, and hence the control term e^f is of order e^. In a very 
similar way, if we assume that the Hamiltonian is three times different iable in (A, 0), 
then the derivative of the control term is also bounded and of order e^. 
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